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Abstract 

Quasi-Monte Carlo (QMC) sampling has been developed for integra¬ 
tion over [0, l] s where it has superior accuracy to Monte Carlo (MC) for 
integrands of bounded variation. Scrambled net quadrature gives allows 
replication based error estimation for QMC with at least the same accu¬ 
racy and for smooth enough integrands even better accuracy than plain 
QMC. Integration over triangles, spheres, disks and Cartesian products of 
such spaces is more difficult for QMC because the induced integrand on a 
unit cube may fail to have the desired regularity. In this paper, we present 
a construction of point sets for numerical integration over Cartesian prod¬ 
ucts of s spaces of dimension d, with triangles (d = 2) being of special 
interest. The point sets are transformations of randomized (t, m, s)-nets 
using recursive geometric partitions. The resulting integral estimates are 
unbiased and their variance is o(l/n) for any integrand in L 2 of the prod¬ 
uct space. Under smoothness assumptions on the integrand, our random¬ 
ized QMC algorithm has variance 0(n~ 1 ~ 2 ' d (\ogn) s ~ 1 ), for integration 
over s-fold Cartesian products of d-dimensional domains, compared to 
0(n -1 ) for ordinary Monte Carlo. 

Keywords: Discrepancy, Multiresolution, Rendering, Scrambled net, Quasi- 
Monte Carlo. 


1 Introduction 

Quasi-Monte Carlo (QMC) sampling is designed for problems of integration over 
the unit cube [0, l] s . Sampling over more complicated regions, such as the tri¬ 
angle or the sphere is more challenging. Measure preserving mappings from the 
unit cube to those spaces work very well for plain Monte Carlo. Unfortunately, 
the composition of the integrand with such a mapping may fail to have even the 
mild smoothness properties that QMC exploits. 

In this paper, we consider quasi-Monte Carlo integration over product spaces 
of the form X s where X is a bounded set of dimension d. We are especially 
interested in cases with d = 2 such as triangles, spherical triangles, spheres, 
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hemispheres and disks. Integration over such sets is important in graphical 
rendering (Arvo et al. 2001). For instance, when A is a triangle, an integral 
of the form”J^7"/( ;:(; i, * 2 ) da^i d ^2 describes the potential for light to leave one 
triangle and reach another. The function / incorporates the shapes and relative 
positions of these triangles as well as whatever lies between them. 


Recent work by Basil and Owen (2014) develops two QMC methods for use 


in the triangle. One is a lattice like construction that was the first construction 
to attain discrepancy 0(log(n)/n) in that space. The other is a generalization 
of the van der Corput sequence that makes a recursive partition of the triangle. 

In this paper, we generalize that van der Corput construction from the unit 
triangle to some other sets. We also replace the van der Corput sequence by 
digital nets in dimension s, to obtain QMC points in X s . The attraction of 
digital nets is that they can be randomized in order to estimate our quadrature 
error through independent replication of the estimate. Those randomizations 
have the further advantage of reducing the error by about 0(n~ 1 ^ 2 ) compared 
to unrandomized QMC, when the integrand is smooth enough. For a survey of 
randomized QMC (RQMC) in general, see L’Ecuyer an d Lemieux| ( [2002 ). For 
an outline of QMC for computer graphics, see Keller (2013). 

We study QMC and RQMC estimates of 
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vol(A) s 

Our estimates are equal weight rules 


fix) Ax. 


1 X s 
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n 


A = - f(xi), where a;» = <j>{ui) 
n z ' 

2=1 


(1) 


for random points Ui £ [0, l] s . The transformation </> maps [0,1) into X and is 
applied componentwise. We assume throughout that vol(A) = 1 whenever we 
are integrating over X , which simplifies several expressions. 

For any / £ L 2 {X S ) we find that Var(A) = o(l/n), so it is asymptotically 
superior to plain Monte Carlo. We also find that for each finite n, scrambled 
nets have a variance bounded by a finite multiple of the Monte Carlo variance, 
uniformly over all / £ L 2 {X S ). 

Our main result is that under smoothness conditions on / and a sphericity 
constraint on the partitioning of X we are able to show that the estimate ([!]) 
attains 

v«(A) = o(fl^d) (2) 

when Ui are certain scrambled digital nets. This variance rate is obtained via 
a functional ANOVA decomposition of the integrand. The case d = 1 in § 
corresponds to the rate for scrambled net RQMC from Owen (1997b). The 


primary technical challenge in lifting that result from d = 1 to d > 1 is to show 
that the composition / o is a well-behaved integrand. 

The statements above are for integration over X s but our proof of the main 
result is for integration over n-= .j X^ where X^ C are potentially different 
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sets of dimension d. Some of our results allow different dimensions dj for the 
spaces X^\ 

The rest of the paper is organized as follows. In Section [2j we give back¬ 
ground material on digital nets and their scrambling. In Section |3j we present 
recursive geometric splits of a region X C and geometric van der Corput se¬ 
quences based on them. Section [4] generalizes those constructions to Cartesian 
products of s ^ 1 such sets. Section [5] presents the ANOVA and multiresolu¬ 
tion analysis of the Cartesian product domains we study. Those domains are 
not rectangular and we embed them in rectangular domains and extend the 
integrands to rectangular domains as described in Section [6j We use both a 
Whitney extension and a Sobol’ extension and give new results for the latter. 
The proof of our main result is in Section [7] Section [8] compares the results we 
obtain to plain QMC and scrambled nets over sd dimensions. 

We conclude this section by citing some related work on QMC over tensor 
product spaces. Tractability results have been obtained for integration over the 
s-fold product of the hypersphere S d = {x £ M' ,+1 | x T x = 1} by Kuo and 
Sloan (2005). |Basu| ( 2014| ) obtained such results for the s-fold product of the 
simplex T d = {x £ [0, l] d | x .i ^ !}■ Those results are non-constructive. For 
s-fold tensor products of S 2 there is a component-by-component construction 
by Hesse et ah (12007). 


2 Background on QMC and RQMC 


Both QMC and ordinary Monte Carlo (MC) correspond to the case with d = 1 
and X = [0,1], Plain Monte Carlo sampling of [0, l] s takes x, ~ U[0, l] s . The 
law of large numbers gives ft —> /i with probability one when / £ L 1 . If also 
f £ L 2 then the root mean square error is a / yfn where a 2 is the variance of 
f(x) for x ~ U[0, l] s . 

QMC sampling improves upon MC by taking a % more uniformly distributed 
in [0, l] s than random points usually are. Uniformity is measured via discrep¬ 
ancy. The local discrepancy of aq,..., x n £ [0, l] s at point a £ [0, l] s is 

1 n 

6(a) = S(a; x lt ..., x n ) = - V l Xi e[n.a) - vol([0, a)). 

n “j' 

The star discrepancy of those points is 

.,*„)= D* n = sup |<5(a)|. 

aG[0,l] d 

The Koksma-Hlawka inequality is 




where Vhk is the s-dimensional total variation of / in the sense of Hardy and 
Krause. For a detailed account of Vhk see Owen ( |2005). Numerous c onstruc- 
tions are known for which D* = O ((log «) s_i /n) (Niederreiter, 1992) and so 
QMC is asymptotically much more accurate than MC when Vhk(/) < oo. 
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2.1 Digital nets and sequences 


Of special interest here are QMC constructions known as digital nets (Niederre- 


iter 


1987 Dick and Pillichshammer 2010). We describe them through a series 


of definitions. Throughout these definitions b ^ 2 is an integer base, s ^ 1 is an 
integer dimension and Z& = { 0 , 1 ,..., 5 — 1 }. 


Definition 1. For kj G N and Cj G Z& for j = 1, , s , the set 



i 


is a 6 -adic box of dimension s. 

Definition 2. For integers m ^ t ^ 0, the points xi,. .. , x^m G [0, l] s are 
a (f, to, s)-net in base b if every 6 -adic box of dimension s with volume b t ~ m 
contains precisely 6 * of the Xi . 

The nets have good equidistribution (low discrepancy) because boxes [0, a] 
can be efficiently approximated by unions of &-adic boxes. Digital nets can 
attain a discrepancy of 0 ((log(n)) s_ 1 /n). 

Definition 3. For integer t ^ 0, the infinite sequence Xi,X 2 , ■ • • G [0, l] s is a 
(f , s)-sequence in base b if the subsequence Xi +r i ,™,..., x^+i^m is a (f , to , s)-net 
in base b for all integers r ^ 0 and to ^ t. 

The (f, s)-sequences (called digital sequences) are extensible versions of (t, to, s)- 
nets. They attain a discrepancy of 0((log(n)) s /n). ft improves to 0((log(n)) s_ 1 /n) 
along the subsequence n = A b m for integers m ^ 0 and 1 ^ A < b. 


2.2 Scrambling 


Here we consider scrambling of digital nets and give several theorems for [0, l) s 
that we generalize to X s . Let a G [0,1) have base b expansion a = J2V= l a kb~ k 
where a*, G ZIf a has two base b expansions, we take the one with a tail of 
0s, not a tail of b — Is. We apply random permutations to the digits ak yielding 
Xk G Zf, and deliver x = x kb~ k ■ There are many different ways to choose 


from Owen (1995). 


the permutations (Owen 2003). Here we present the nested uniform scramble 


In a nested uniform scramble, X\ = 7r.(ai) where 7r. is a uniform random 
permutation (all b\ permutations equally probable). Then x 2 = 7 r. ai (a 2 ), £3 = 
’^.a 1 a 2 ( a 2 ) and Xk +1 = TT. ai a 2 ...a k (afc+i) where all of these permutations are 
independent and uniform. Notice that the permutation applied to digit dk+i 
depends on the previous digits. A nested uniform scramble of a = (a \,..., a s ) G 
[ 0 , l) s applies independent nested uniform scrambles to all s components of a, so 
that Xj tk+ 1 = ■Kj.a :il aj 2 ,...,a jk {aj,k+\)- A nested uniform scramble of a\,. .., a n G 
[ 0 , l) s applies the same set of permutations to the digits of all n of those points. 
Propositions [l] and [ 2 ] are from Owen (1995). 
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Proposition 1. Let a £ [0, l) s and let x be the result of a nested uniform 
random scramble of a. Then x ~ U[0, l) s . 

Proposition 2. If the sequence a 1; ..., a n is a (t,m, s)-net in base b, and x^ 
are a nested uniform scramble of at, then Xi are a (t,m, s)-net in base b with 
probability 1. Similarly if is a (f, s)-sequence in base b, then Xi is a (t,s)~ 
sequence in base b with probability 1. 

In scrambled net quadrature we estimate p = Jj 0 ^ f{x) dx by 

1 ” 

£= (3) 

2=1 


where Xi are a nested uniform scramble of a digital net a^. 

It follows from Proposition 0 that E(/i) = p for / £ L^O, l) s . When 
/ £ L 2 [0, l) s we can use independent random replications of the scrambled 
nets to estimate the variance of p. If Vhk(/) < oo then we obtain Var(/i) = 
0(log(n) 2 ( s-1 )/?i 2 ) = 0(n~ 2+e ) for any e > 0 directly from the Koksma-Hlawka 
inequality. Surprisingly, scrambling the net has the potential to improve accu¬ 
racy: 

Theorem 1. Let f : [0, l] s —► R with continuous dxi d .. dx f ■ Suppose that Xi 
are a nested uniform scramble of the first n = A b m points of a (t, s)-sequence in 
base b, for A £ {1, 2,..., b — 1}. Then for p given by §. 

Var W = 0( l0g( " r ‘ )=0("- 3+ ‘) 


as n —)■ oo for any e > 0. 


Proof. Owen (1997b) has this under a Lipschitz condition. Owen (2008) removes 
that condition and corrects a Lemma from the first paper. □ 


Smoothness is not necessary for scrambled nets to attain a better rate than 
Monte Carlo. Bounded variation is not even necessary: 

Theorem 2. Let xi, ..., x n be a nested uniform scramble of a (t, to, s)-net in 
base b. Let f £ L 2 ([0, l] s ). Then for p given by ([3|, 

v arW = o(t) 

os n —> oo. 


Proof. This follows from Owen (1998). The case t = 0 is in Owen (1997a). □ 


The factor log(ro) s_1 is not necessarily small compared to n 3 for reasonable 
sizes of n and large s. Informally speaking those powers cannot take effect for 
scrambled nets until after they are too small to make the result much worse 
than plain Monte Carlo: 


5 
















Theorem 3. Let aq,..., x n be a nested uniform scramble of a (t, to, s)-net in 
base b. Let f £ L 2 ([ 0, l] s ) with Var (f(x)) = a 2 when x ~ U[0, l] s . Then for (i 
given by ([3|, 

If t = 0, then Var(/t) ^ ea 2 /n = 2.718<r 2 /n. 


Proof The first result is in Owen (1998), the second is in Owen (1997a). □ 


3 Splits and geometric van der Corput sequences 


The van der Corput sequence is constructed as follows. We begin with an 
integer i ^ 0. We write it as i = J2kLi a k(i)b k ~ 1 for digits a^i) £ Z b- Define 
the radical inverse function <j>b(i) = a k{i)b~ k £ [0,1]. The van der Corput 

sequence in base b is Xi = <j>b(i — 1) for i ^ 1. It is a (0, l)-sequence in base b. 
The original sequence of van der Corput. (1935a] b) was for base 6 = 2. Any n 
consecutive van der Corput points have a discrepancy 0(\og{n)/n) where the 
implied constant can depend on b. 

The lowest order base 6 digit of i determines which of 6 subintervals [a/6, (a+ 
1) /b) will contain aThe second digit places Xi into one of 6 sub-subintervals of 
the subinterval that the first digit placed it in, and so on. Basu and Owen (2014| 
used a base 4 recursive partitioning of the triangle to generate a triangular 
van der Corput sequence. Discrepancy in the triangle is measured through 
equidistribution over trapezoidal subsets (Brandolini et ah, 2013). Triangular 


van der Corput points have trapezoidal discrepancy of 0(n i/2 ). 


3.1 Splits and recursive splittings 

We begin with a notion of splitting sets. Splits are like partitions, except that 
we don’t require empty intersections among their parts. 

Definition 4. Let X C have finite and positive volume. A 6-fold split of X is 
a collection of Borel sets X a for a £ Zb with X = U„I.g X a , vol(A a ) = vol(X)/b 
for a £ Zb, and vol(T a n X a >) = 0 for 0 ^ a < a' < b. 

In all cases of interest to us, any overlap between X a and X a / for a a' 
takes place on the boundaries of those sets. The unit interval [0,1) is custom¬ 
arily partitioned into subintervals [a/6, (a + l)/6) in QMC. Handling X = [0,1] 
requires awkward exceptions where the rightmost interval is closed and all oth¬ 
ers are half-open. For general closed sets X it could be burdensome to keep 
track of which subsets had which parts of their boundaries. Using splits allows 
one for example to divide a closed triangle into four congruent closed triangles. 
Of course a partition is also a valid split. Our preferred approach uses a ran¬ 
domization under which there is probability zero of any sample point appearing 
on a split boundary. 
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6 = 2 


6=3 


6=4 



Figure 1: Splits of a triangle A” for bases 6 = 2, 3 and 4. The subtriangles Xj 
are labeled by the digit j £ Zj,. 

Figure |T| shows a triangle X = ABC split into subtriangles. The left panel 
has 6 = 2 subtriangles, the middle panel has 6 = 3 and the right panel has 6 = 4. 
For 6 = 2, the vertex labeled ‘A’ is connected to the midpoint of the opposite 
side. The subset ABCo is the one containing ‘B’. In each case, the new ‘A’ is 

the mean of the old ‘B’ and ‘C’. The new ‘B’ is to the right as one looks from 

the new A towards the center of the triangle, and the new ‘C’ is on the left. An 
algebraic description is more precise: Using lower case abc to describe the new 
ABC , the case described above (base 2 and digit 0) has 

a\ /0 1/2 l/2\ 

6 = 1 0 0 

c) \0 1 0 / 

Similar rules apply to the other bases. From such rules we may obtain the 
vertices of a split at level k by multiplying the original vertices by a sequence 
of k 3 x 3 matrices operating on points in the plane. 

Definition 5. Let X C have finite and positive volume. A recursive 6-fold 
split of A is a collection X of sets consisting of X and exactly one 6-fold split of 
every set in the collection. The members of X are called cells. 

The original set X is said to be at level 0 of the recursive split. The cells 
Xo,..., Xb-i of X are at level 1. A member of a recursive split of X is at 
level k ^ 1 if it arises after k splits of X. The cell X ai: a 2 is the subset of X ai 
corresponding to a = 02 and similarly, an arbitrary cell at level k ^ 2 is written 
Xai ,0.2 ,... for a.j £ Z ib . 

We will need to enumerate all of the cells in a split X. For this we write 
t = Y^ k j =1 a j^~ 1 £ 2 b k and then take X^t) = X aii a 2 ^ ^ ak The cells in the split 
are now X^ t \ for k £ N and t £ Z b k , with A( 0 ,o) = A. 

Figure [2] shows the first few levels of recursive splits for each of the splits from 
Figure]!] The base 3 version has elements that become arbitrarily elongated as k 
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2 6 Decomposition 


3 3 Decomposition 


4 3 Decomposition 



Figure 2: The base b splits from Figure [l] carried out to k = 6 or 3 or 4 levels. 


increases. That is not desirable and our best results do not apply to such splits. 
The base 4 version was used by|Basu and Owen (2014) to define a triangular van 
der Corput sequence. The base 2 version at 6 levels has a superficially similar 
appearance to the base 4 version at 3 levels. But only the latter has subtriangles 
similar to X. A linear transformation to make the parent ABC an equilateral 
triangle yields congruent equilateral subtriangles for b = 4, while for b = 2 one 
gets isoceles triangles of two different shapes, some of which are elongated. 


3.2 Splitting the disk and spherical triangle 

The triangular sets were split in the same way at each level. Splits can be more 
general than that, as we illustrate by splitting a disk. A convenient way to define 
a subset of a disk is in polar coordinates via upper and lower limits on the radius 
and an interval of angles (which could wrap around 27r = 0). If one alternately 
splits on angles and radii, the result is a decomposition of the disk into cells 
of which some have very bad aspect ratios, especially those near the center. 


Beckers and Beckers (2012) define the aspect ratio of a cell as the ratio of the 


length of a circular arc through its centroid to the length of a radial line through 
that centroid. They show decompositions of the disk into b cells with aspect 
ratios near one for b as large as several hundred. But their decompositions are 
not recursive. Figure [3] shows eight levels in a recursive binary split of the disk 
into cells. A cell with aspect ratio larger than one is split along the radial line 
through its centroid. Other cells are split into equal areas by an arc through 
the centroid. 

A spherical triangle is a subset of the sphere in K 3 bounded by 3 great circles. 
By convention, one only considers spherical triangles for which all internal angles 
are less than 7r radians. The spherical triangle can be split into b = 4 cells like 
the rightmost panel in Figure [I] If one does so naively, via great circles through 
midpoints of the sides of the original triangle, the four cells need not have equal 
areas. Song et al. (2002]) present a four-fold equal area recursive splitting for 






















































Adaptive binary splits of a disk 







Figure 3: A recursive binary equal area splitting of the unit disk, keeping the 
aspect ratio close to unity. 


spherical triangles, but their inner boundary arcs are in general small circles. If 
one wants a recursive splitting of great circles into great circles, then it can be 
done by generalizing the 6 = 2 construction of the leftmost panel in Figure [l] A 
great circle with vertices ABC can be split into two by finding a point P on BC 
so that APB has half the area of ABC, using the first step of Arvo’s algorithm 


Arvo (19951. 


3.3 Geometric van der Corput sequences 

Given a set X and a recursive splitting of it in base 6 we can construct a 
geometric van der Corput sequence for X . The integer i is written in base b as 
i = Y^kLi a k{i)b k ~ 1 - To this i we define a sequence of sets 


Xj:K {i) ‘ 

Then x t is any point in n^ =1 Aj./f. The volume of Xi K is b~ K which converges 
to 0 as K —> oo. For most of the constructions we are interested in, each a % is a 
uniquely determined point. For decompositions with bad aspect ratios, like the 
base 3 decomposition in Figure [2j some of the set sequences converge to a line 
segment. For instance if i £ {0,1, 2}, then a*.(i) = 0 for k ^ 1 and that infinite 
tail of zeros leads to a point Xi on one of the edges of the triangle. 

To get a unique limit Xi, we use the notion of a sequence of sets converging 
nicely to a point. Here is the version from Stromberg (1994). 
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Definition 6. The sequence Sk £ of Borel sets for k £ N converges nicely 
to x £ R. d as k —» oo if there exists a < oo and c?-dimensional cubes Ck such 
that x £ Ck , Sk C Cfc, 0 < vol(Cfc) ^ avol(<Sfc), and lim^oo diam(<Sfc) = 0. 

A sequence of sets that converges nicely to x cannot also converge nicely to 
any x' ^ x. We generally assume the following condition. 

Definition 7. A recursive split X in base b is convergent if for every infinite 
sequence aq, oq, 03 , • • • € Zf,, the cells A Qli02) ...« x converges nicely to a point as 
K —> 00 . That point is denoted limA'->oo X aija2t ^ taK . 

In a geometric van der Corput sequence, we take a convergent recursive split 
and choose 

Xi = A ’“i(i-l).“2(i-l),...,OK(i-l),p,0, ... ,0 

V" y 

M 

where K is the last nonzero digit in the expansion of i — 1 and there are M ^ 1 
zeros above. For the base 4 triangular splits, aq is simply the center point of 
‘^ai(i-i),a 2 (*-l),...,ajr(i-i)- For b = 2, Xi is an interior but noncentral point. The 
recursive split for b = 3 is not convergent. 

Definition 8. Let X be a recursive split of X € in base b. Then X satisfies 
the sphericity condition if there exists C < oo such that diam(A’ aii ... jafc ) ^ 
Cb~ k ! d holds for all cells Af 0l) ... )0fc in X. 

A recursive split that satisfies the sphericity condition is necessarily conver¬ 
gent. The constant C can be as low as 1 when d = 1 and the cells are intervals. 
The smallest possible C is usually greater than 1. We will assume without loss 
of generality that 1 ^ C < oo. 

Definition 9. Given a set X C and a convergent recursive split X of A 
in base b , the X-transformation of [0,1) is the function <j> = cf>x : [0,1) — > X 
given by 4>(x) = lim k^-oo X XltX2 ^.. tXK where x has the base b representation 
O.X 1 X 2 ■ ■ • • If 2 : has two representations, then the one with trailing Os is used. 


4 Geometric nets and scrambled geometric nets 

Let A be a bounded subset of with finite nonzero volume. Here we define 
digital geometric nets in X s via splittings. It is convenient at this point to 
generalize to an s-fold Cartesian product of potentially different spaces, even 
though they may all be copies of the same A. 

For s € N, we represent the set {1,2,..., s} by l:s. For j £ l:s we have 
bounded sets A^ C with vol(A^^) = 1. For sets of indices u C l:s, the 
complement l:s — u is denoted by —u. We use m| for the cardinality of u. The 
Cartesian product of A^- 1 for j £ u is denoted X u . A vector x £ X 1:s has 
components Xj £ X^k The vector in A“ with components x 3 for j £ u is 
denoted x u . 
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A point in X 1:s has ]Cj=i dj components. We write it as x = (xi,x 2 ,..., x s ). 
The components in this vector of vectors are those of the Xj concatenated. The 
notation Xj for dj consecutive components of x is the same as we use for the 
j’th point in a quadrature rule. The usages are different enough that context 
will make it clear which is intended. 

Definition 10. For j = 1,... ,s, let Xj be a recursive split of X^> in a common 
base b. Denote the cells of Xj by Xjj kt j for k £ N and t £ Z b k. Then a 6-adic 
cell for these splits is a Cartesian product of the form n;=i *j,( M .) for integers 
kj ^ 0 and tj € Z^. 

Definition 11. Let X^) c M dj have volume 1 for j £ l:s and let Xj be a 
recursive split of XA) in a common base b. For integers m ^ t ^ 0, the points 
xi,... ,Xbm £ X 1:s are a geometric {t,m,s )-net in base b if every 6-adic cell 
of volume contains precisely fe* of the »$. They are a weak geometric 

( t , m, s)-net in base b if every 6-adic cell of volume b t ~ m contains at least b f of 
the Xi. 

Some of the 6-adic cells can get more than b l points of a weak geometric 
(t, m , s)-net because the boundaries of those cells are permitted to overlap. 

Proposition 3. Let a i,..., a n be a ( t, m, s)-net in base b. Let u±,... ,u n be a 
nested uniform scramble of a±,... ,a n . For j £ l:s, let Xj be a recursive base b 
split of the unit volume set X^ C with transformation (j>j. Then Zi = <j>(a,i ) 
(componentwise) is a weak geometric (t,m, s)-net in base b and Xi = (f(uf) 
(componentwise) is a geometric (t,m, s)-net in base b with probability one. 

Proof. In both cases the transformation applied to half open intervals places 
enough points in each 6-adic cell to make those points a weak geometric (f, m, s)- 
net. The result for scrambled nets follows because each x, is uniformly dis¬ 
tributed and 

vol (Aj i(M) n^,(M')) =0 

for all j £ l:s, k £ Z b k and 0 ^t<t'<b k . □ 

4.1 Measure preservation 

Let <j) : [0,1) —> X C be the function that takes a point in the unit interval 
and maps it to X according to the convergent recursive split X. We show here 
that <)> preserves the uniform distribution. This is the only section in which we 
need to distinguish Lebesgue measures of different dimensions. To that end, we 
use Ai for Lebesgue measure in R and A d for R d . 

Proposition 4. Let X C R d with vol(A) = 1. Let X be a convergent recursive 
split of X in base b ^ 2. Let <f) be the X -transformation of [0,1) and let AC X 
be a Borel set. Then 

where (j)~ x {A) = {x £ [0,1) | cj>(x) £ A}. 
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Proof. First, suppose that A = X ail a 2r .., ak for aj G Z b . Then </> X (A) = 
[: t/b k , (t + 1 )/b k ) for some t G Z b fc, and so 


A d (A) 


^- Al l 


vLfofc’ b k )) 


Ai (r 1 ^)). 


Now let A be any Borel subset of X. Given e > 0, there exists a level k\ ^ 1 
and n cells Bi at level k\ of X such that A C ULi A and A d (Ur =1 B< \ Ai) < c- 
Similarly, there exists a level k 2^1 and to cells C,; at level &2 of X such that 
UHi CiC A and A d (A \ U'=i C*) < e. Thus we get, 


A d (A) = X d (^ 0 Bi) - A d ^ U B * \ A ) > Ad( 0 ^ e 

= £ Ad(Bi) - c = Ai(0 _1 (Bi)) ~ e = A, (V 1 ( U B, 

i = l i=l A Vi—1 

> Ai(<^ 1 (-4)) - e, 


where the second equality follows because (f> is bijective and therefore <f~ 1 (A ) C 
0 _ 1 (Ur=i Bi). Similarly, X d (A) < Ai(0 _ 1 (A)) + e. Since e was arbitrary we 
have the proof. □ 

Measure preservation extends to the multidimensional case. The next propo¬ 
sition combines that with uniformity under scrambling. 

Proposition 5. Let X^ C R dj with vol(dfo^) = 1 for j G l:s have convergent 
recursive splits Xj in bases bj ^ 2 with corresponding transformations <fj. Let 
a G [0, l) s and let x„ be a base 6., nested uniform scramble of a.;. Then d>(x] = 

Proof. By Proposition [I] x is uniformly distributed on [0, l) s . By Proposition]^] 
(f> preserves uniform measure. Thus <pj(xj ) are independent XJ(X < A'} random 
elements. □ 


4.2 Results in L 2 not requiring smoothness 

Some of the basic properties of scrambled nets go through for geometric scram¬ 
bled nets, without requiring any smoothness of the integrand. They don’t even 
require that the same base be used to define both the transformations and the 
digital net. 

Theorem 4. Let XA) c with vol(Afh)) = 1 for j = 1,..., s. Let Xj be 
a convergent recursive split of X^ in base bj ^ 2 with transformation 4>j. Let 
Hi,..., u n be a nested uniform scramble of a (f, to, s)-net in base b ^ 2 and let 
Xi = <f)(ui) componentwise. Then for any f G L 2 (X 1:s ), 

VarM=o(i) 

as n —> oo. 
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Proof. Since / £ L 2 (X 1:s ) we have / o </> £ L 2 [ 0, l] 5 . Then Theorem [ 2 ] applies. 

□ 


Theorem 5. Under the conditions of Theorem [7} 

where a 2 = Var (f(x)) for x ~ U(A 1:S ). If t = 0, then Var(/i) ^ ea 2 /n = 
2.718 a 2 /n. 

Proof. Once again, / o <f> £ L 2 [0, l] s . Therefore Theorem [ 3 ] applies. □ 


5 ANOVA and multiresolution for X 1:s 


There is a well known analysis of variance (ANOVA) for [0, l) s . Here we present 
the corresponding ANOVA for X 1 ' 8 . Then we give a multiresolution of L 2 (X 1:s ) 
adapting the base b wavelet multiresolution in Owen (1997a) for [0,1). 


5.1 ANOVA of X 1:s 

For / £ L 2 (X 1:s ) the ANOVA decomposition provides a term for each u C l:s. 
These are defined recursively via 



The function f u represents the ‘effect’ of x 3 for j £ u above and beyond what 
can be explained by lower order effects of strict subsets v C u. While f u 
is a function defined on X 1:s its value only depends on x u . By convention 
fz{x) = J x i:s f(x)dx = 11 for all x. We define variances a 2 = J x i-„ f u (x) 2 dx 
for |u| > 0 and cr| = 0. The ANOVA decomposition satisfies X]|«|>o cr u = 0-2 
where a 2 = f x i , s {f(x) — n) 2 dx. It also satisfies f(x) = X)uCi- s fu{ x ) by the 
definition of fi :s , wherein a 0-fold integral of a function leaves it unchanged. 


5.2 Multiresolution 

We begin with a version of base b Haar wavelets adapted to X C using a 
recursive split X of A in base b ^ 2. Recall that the cells at level k of a split 
are represented by one of X^ k t ) for 0 ^ t < b k . Those cells are in turn split at 
level k + 1 via 

b -1 

k,t t c)j wheie ^(k,t,c) ‘^(fc+i,6t-i-c)* 

c=0 
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The multiresolution of X in terms of X has a function <p{x) = 1 for all x £ X 
as well as functions 


i’ktc = 6 (fc+1)/2 l*e* (M , o) - 6 (fe - 1)/2 


(5) 


where 7Vfc tc and ITfct are indicator functions of the given narrow and wide cells 
respectively. The scaling in ([5]) makes the norm of ipktc independent of k : 
f ipl tc (x) dx = (b — 1 )/b. 

For /i,/ 2 £ L 2 (X) define the inner product {fi, f%) = f x fi( x )f 2 (x) da;. 
Then let 

K b k -1 6-1 

Jk(x) = (f,<p)<p(x)+J2 Z J2(f^ktc)i>ktc(x). 

k—1 t—0 c—0 

For x belonging to only one cell at level K + 1, as almost all x do, Jk(x) 
is the average of / over that cell. By Lebesgue’s differentiation theorem, local 
averages over sets that converge nicely to x satisfy 


lim 

K—t oo 


Isk / W d:E 

vol (5if) 


/(*), a.e. 


for / £ L 1 (M. d ). So if X is convergent, then lim/f-^oo /if (x) = f(x) almost 
everywhere. Thus we may use the representation 


oo 6 fc —1 6—1 

/(*) = (/, ¥>M*) + EEE^’ kktc)kktc(x). (6) 

k=1 t =0 c—0 

Equation Q resembles a Fourier analysis with basis functions tp and ipktc- Un¬ 
like the Fourier case, the functions ipktc and 'P’ktc' are not orthogonal. Indeed 
J2ce Zb ^ fetc = 0 a - e " Non-orthogonal bases that nonetheless obey © are known 
as tight frames. 

We may extend © to the multidimensional setting by taking tensor prod¬ 
ucts. For j £ l:s, let X 1 ^ C have recursive split X ? - in base b ^ 2. Let the 
basis functions be ipj and ipj(ktc) with narrow and wide cell indicators Njktc and 
Wjkt • For u C l:s, let k £ kjUl have elements kj ^ 0 for j £ u. Similarly let r 
have elements tj £ Z^ and 7 have elements Cj £ Z&, both for j £ u. Then for 
x £ X 1:s define 


' l Pu.KTy{x) — ~W' l Pjkjt : jCj(Xj) E[v 5 j'( a: j)- (7) 

j£u jgu 

Our multiresolution of L 2 (X 1:S ) is 


/(*) = z z z Z^ UKT^l fy'lfiuKT'Y (®) 

itCl:s k|u r|n,Ac 7 |ii 

= ^ZZZZ(^ UKT'yi f')1pUKT'y(x'). 

|-zx|>0 K.|xt 7|iz 

The sum over k is over all possible values of k given the subset it. The other 
sums are similarly over their entire ranges given the other named variables. 
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5.3 Variance and gain coefficients 

Here we study the variance of averages over scrambled geometric nets. We start 
with arbitrary points a.j £ [0, l) s . For now, they need not be from a digital net. 
They are given a nested uniform scramble, yielding points Ui £ [0, l) s . Those 
points are then mapped to x % £ X 1:s using recursive splits in base b. 

It follows from Proposition [5] that 

✓ 1 n n 

Var (fi) = E feEEEEEEE E E 

2=1 i' — 1 |ti|>0 K,\u t\u,K 7| U I tfc 7 I >0 k'\u' t'\u' ,K f 7 7 | u' 

(/> 'IpUKT’y) (/? < 0W / « / T / 7 / )V J UK.T'y (*Ti) 1^u , K r T lr y l ) ]■ 


This formula simplifies due to properties of the randomization. Lemma 4 from 


Owen (1997a) shows that if u 7 ^ u’ or k 7 ^ k! or r 7 ^ r', then, 


^‘(ipUKT'y ) 'Ipu'(*^i' ) ) — 0- 


( 8 ) 


Consequently 

Var(/t) = ^ E Var (-E^^)J, 

|n|>0 /t|n ^ *=1 ^ 

where 

(«^) = ££</ J *0 UKT 7 )^ UKT'y (®) 
r|ii,/c 7|it 

with z/ 0i () = /r. The function z/ UK is constant within elementary regions of the 
form 

n n 

jeu jgu 

for 0 < tj < b kj and 0 ^ Cj < b. 

Define 

o-L = / «2 B (*)d®. 

Jx 1 -* 

The multiresolution-based ANOVA decomposition is 



(9) 


which follows from the orthogonality in ([ 8 ]). 

The equidistribution properties of ai,..., a n determine the contribution of 
each v UK to Var(/t). Write a* = (an ,..., aj S ) and define 

= ^ 1 [b fc + 1 a ij 'J = L bfc + la i'.iJ _ ^ L bfca ij J = L bfca i'jJ ) ' 
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For each |u| > 0 and n £ Nl“l define 


r “’ K ~ „ E E II 


kn * 


i—1 i' — l jEv 


It follows from Theorem 2 of Owen (1997a I that 


1 


Var(A) = - E y~! r “.^ c 


|'Ll,| >0 K,\u 

We must have T UiK ^ 0 because Var(/i) ^ 0. In plain Monte Carlo sampling, 
Var(/t) = a 2 /n which corresponds to all T UjK = 1 (compare ([9])). The T U>K are 
called ‘gain coefficients’ because they describe variance relative to plain Monte 
Carlo. If the points eq are carefully chosen, then many of those coefficients can 
be reduced and an improvement over plain Monte Carlo can be obtained. 

If ai,...,a n are a (t, m, s)-net in base 6, we can put bounds on the gain 
coefficients using lemmas from Owen (1998). In particular, F UjK = 0 if m — t ^ 

|it | + |ac|, and otherwise 

'6 + 1 ' ' 

. 6 - 1 . 


r„, K < 6* 


Thus finally we have, 


I/' 

Var (fi) < — 
n 


6+1 

6-1 


E E < 

|w|>0 |«| + |n|>m —t 


( 10 ) 


Equation (101 shows that the scrambled net variance depends on the rate 
at which a 2 K decay as |/t| + |u| increases. For smooth functions on [0, l) s they 
decay rapidly enough to give Var(/i) = 0(log(n) s_1 /n 3 ). To get a variance rate 
on X 1:s we study the effects of smoothness on a 2 K for d-dimensional spaces X. 


6 Smoothness and Extension 

Our main results for scrambling geometric nets require some smoothness of 
the integrand. We also use some extensions of the integrand and its ANOVA 
components to rectangular domains. 

Let / be a real-valued function on X C R m . The dimension m will usually 
be d x s, for an s-fold tensor product of a d-dimensional region. For v C l:m, 
the mixed partial derivative of / taken once with respect to Xj for each j £ v 
is denoted d v f. By convention d 0 f = /, as differentiating a function 0 times 
leaves it unchanged. 

6.1 Sobol’ extension 

We present the Sobol’ extension through a series of definitions. 
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Definition 12. Let X C R m for m G N. The function / : X — > R. is said to be 
smooth if d 1:m f is continuous on X. 

Definition 13. Let X C R m . The rectangular hull of X is the Cartesian 
product 

m 

rect(T’) = ]^J[inf{a;j | x G A}, sup{#j | x G X}], 

j =i 

which we also call a bounding box. For two points a, b G R m we write rect[a, b] 
as a shorthand for rect[{a, b}}. 

For later use, we note that 

diam(rect(T’)) ^ Vd x diam(T’), (11) 


for lcl d 

Definition 14. A closed set X C R m with non-empty interior is said to be 
Sobol’ extensible if there exists a point c G X such that z G X implies rect[c, z] C 
X. The point c is called the anchor. 

Figure [4] shows some Sobol’ extensible regions. Figure [5] shows some sets 
which are not Sobol’ extensible, because no anchor point exists for them. Sets 
like the right panel of Figure [5] are of interest in QMC for functions with inte¬ 
grate singularities along the diagonal. 

Proposition 6. If X^ C M. dj is Sobol’ extensible with anchor Cj for j = 

1 ,..., s, then n-=i XU) 

is Sobol’ extensible with anchor c= (ci, ..., c s ). 

Proof. Suppose that x G n;= ^ X^\ We write x as (xi, ... , x s ) where each 
Xj G X^\ Then rect(c, x) C rec t(c 7 , x 7 ) G IIj=i ■ $3 

Given points x,y G K m and a set u C l:m, the hybrid point x u :y_ u is the 
point z G R m with Zj = Xj for j G u and Zj = yj for j ^ u. We will also require 
hybrid points x u :y v :z w whose j’th component is that of x or y or z for j in u 
or v or w respectively, where those index sets partition 1 :m. 

A smooth function / can be written as 

/(*) = 51 [ d u f(c_ u :y u )dy u ( 12 ) 


where Jj c x j denotes ± f Tect r c x j ■ The sign is negative if and only if Cj > Xj 
holds for an odd number of indices j G u. The term for u = 0 equals /(c) 
under a natural convention. Equation (121 is a multivariable version of the 
fundamental theorem of calculus. For to = 1 it simplifies to f{x) = f(c ) + 

Sc f(y) d v- 
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Figure 4: Sobol’ extensible regions. At left, X is the triangle with vertices 
(0,0), (0, y/2), (y/2, 0) and the anchor is c = (0,0). At right, A is a circular disk 
centered its anchor c. The dashed lines depict some rectangular hulls joining 
selected points to the anchor. 


Definition 15. Let / be a smooth real-valued function on the Sobol’ extensible 
region X C R m . The Sobol’ extension of / is the function / : R m —> R given by 

f(x) = [ dU f(. c -u-yu) l a_u.y u £X&y u (13) 

where c is the anchor of X. 


The Sobol’ extension can be restricted to any domain X' with X C X' C R m . 
We usually use the Sobol’ extension for / from X to rect(A). This extension 


was usee 

in Sobol’ 

( 

1973 

but not explained there. An account of it appears 

in 

Owen 

(2005 

2006 

). For x £ X the factor l c _ u: j /ii gar is always 1, making 


f{x) = f(x) so that the term “extension” is appropriate. 

Two simple examples serve to illustrate the Sobol’ extension. If X = [0,1] 
and f(x) = x on 0 < x < 1, then the Sobol’ extension of / to [0, oo) is f(x) = 
min(x,l). If X = [0, l] 2 and f{x) = X\X 2 on X, then the Sobol’ extension of 
/ to [0,oo) 2 is f(x) = min(a;i,l) x min(a; 2 ,l). The Sobol’ extension / has a 
continuous mixed partial derivative <9 1:m / for x in the int erior of X and also 


in the interior of R m \ X where d 1:m f = 0 (Owen 2005). As our examples 
show, d u f for |u| > 0 may fail to exist at points of the boundary dX. A 
Sobol’ extensible X C R m has a boundary of m-dimensional measure 0, so 
when forming integrals of d u f we may ignore those points or simply take those 
partial derivatives to be 0 there. 

The Sobol’ extension has a useful property that we need. It satisfies the 
multivariable fundamental theorem of calculus, even though some of its partial 
derivatives may fail to be continuous or even to exist everywhere. We can even 
move the anchor from c to an arbitrary point z. 
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Figure 5: Non-Sobol’ extensible regions. At left, X is an annular region 
centered at the origin. At right, X is the unit square exclusive of an e-wide 
strip centered on the diagonal. 


Theorem 6. Let X C R m be a Sobol’ extensible region and let f have a con¬ 
tinuous mixed partial d 1:m f on X. Let f be the Sobol’ extension of f and let 
z G R m . Then 

/(*)= E I dU f( z -xy u )dy u - (14) 

uCl:m d[z u ,x u ] 

Proof. Define 

g(x) = ^2 [ d v f(z- v :t v )dt v (15) 

uCl ■. m P z '” x d 


where d v f is taken to be 0 on those sets of measure zero where it might not 
exist when |u| > 0. We need to show that g(x) = f(x). Now 

f(z— v 't v ) = I d%f{c_ u :y u )X(c_ u :y u )dy u , (16) 

uCl:m 1 [ c «! 2 wn-u : *iinu] 


where for typographical convenience we have replaced l.^x by X(-). The sub¬ 
script in dy makes it easier to keep track of the variables with respect to which 
that derivative is taken. Substituting (16) into (151, we find that g(x) equals 


m 


vGl :m 1 


E 


iGl:m 


' [c„,z„ n -».:*unu] 


dff(c_ u :y u )X{c_ u :y u ) dy v 


E 

= E / d t E [ d yf( c -x-yu) x ( c -u--yu)dy v 

vCl: L uDv J [ c ^’t u ] 

= E 


d 


d. t v 


E 


U-Vu— IL-Vu — v'^v) dy u _ v d t v . 


vC.l:m ~ 
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Now we introduce w = u — v and rewrite the sum, getting 

E E / [ dy +v f(c_ w _ v :y w :t v )X(c_ w _ v :y w :t v )dy w dt v 

ttCrratC-ed*"'*"] 

= E E / d y [ d™f{c- w - v :y w :t v )X(c- w - v :y w :t v )dy w dt v . 

n,Cl,nvC- w J l^^] 4 «,*»] 

Any term above with v ^ 0 vanishes. Therefore 

9(x) = J2 / dyf{c- w -y w )X{c_ w -.y w )dy w dt e! 

J[o w ,t w ] 


= E 


dyf{c- w :y w )X(c_ w \y w )dy u 


wC.l:m ' 


= /(*)• 


a 


6.2 Whitney extension 


Here we assume that A is a bounded closed set with non-empty interior, not 
necessarily Sobol’ extensible. Sobol’ extensible spaces may fail to have a non¬ 
empty interior, but outside such odd cases, they are a subset of this class. 
Non-Sobol’ extensible regions like those in Figure [5] are included. To handle 
domains X of greater generality, we require greater smoothness of /. 

Let k £ N m be any multi-index with \k\ = k\ + ... + k m < m. We denote 
the fc-th order partial derivative as 


D k f(x) 


dx'l 1 ■ ■ ■ Cfom 


/ Ol, • ■ 


Definition 16. A real-valued function / on A C R m is in C m {X) if all partial 
derivatives of / up to total order m are continuous on X. 

Whitney’s extension of a function in C m (X) to a function in C m (rect(A)) 
is given by the following lemma. 


Lemma 1. Let f £ C m {X) for a bounded closed set X C R m with non-empty 
interior. Then there exists a function f £ C m (rect(A’)) with the following 
properties: 

1. f(x) = f(x) for all x £ X, 

2. D k f{x) = D k f(x) for all |fc| < m and x £ X, and 


3. f is analytic on rect(A) \ X. 

Proof. The extension we need is the one provided by Whitney (1934). A function 
in C m (X) in the ordinary sense is a fortiori in C m (X) according to Whitney’s 
definition. We use the restriction of Whitney’s function to the domain rect(A). 

□ 
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We will need one more condition on X. We require the boundary of X to 
have m-dimensional measure zero. Then Theorem [6] in which the fundamental 
theorem of calculus applies to /, holds also for the Whitney extension. 

6.3 ANOVA components of extensions 

Here we show that the ANOVA components of our smooth extensions are also 
smooth. We suppose that each X' J C M. dj and we let m = The 

Cartesian product X 1:s is now a subset of R m . 

Lemma 2. Let f be a smooth function on Sobol’ extensible X 1:s C R m and for 
u C l:s let f u be the ANOVA component from Q. Then f u is smooth on X 1:s . 

Proof. We prove this by induction on |tt|. Let \u\ = 0, that is u = 0. Then 
f u (x) = J x i:s f(x)dx which is a constant p and is therefore smooth on X 1:s . 
Let us suppose that the hypothesis holds for |u| = k — 1 < s and we shall prove 
it for |u| = k. 

Fix any u C T.s such that |u| = k. By 0 we have, 
fu(x) = f f(x) dx_ u - 

W<ZU 

using the fact that f w (x) does not depend on Xj for j ^ w. Each term in 
the summation is f w for |ui| < k — 1 and is therefore smooth by the induction 
hypothesis. So we only need to show that the first term is smooth. Fix any 
v C 1 :m. Now since / is smooth, d v f{x) is continuous on X 1:s and hence 
applying Leibniz’s integral rule we have, 

d v f f(x)dx_ u = [ d v f(x)dx_ u . 

Now the right hand side is the integral of a continuous function and is therefore a 
continuous function. Thus the induction hypothesis hold for |u| = k completing 
the proof. □ 

Lemma 3. Let f £ C m (X 1:s ) for a bounded closed set X 1:s £ and for 
u C l:s let f u be the ANOVA component in Q. Then f u £ C m (X 1 ' s ). 

Proof. The proof goes along the same lines as Lemma [2j We replace v in that 
proof by any multi-index L with \£.\ < m. Now since / is smooth Dg,f(x) is 
continuous on X 1:s and hence applying the Leibniz’s integral rule we have, 

De f f{x)dx_ u = [ D t f(x)dx_ u . 

JX-v- Jx -« 

Now the right hand side is the integral of a continuous function over certain 
variables and is therefore a continuous function. □ 
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Now for a smooth function / defined on a product X 1:s of Sobol’ extensible 
sets, or on a product of more general spaces but with the smoothness required 
for a Whitney extension, there exists an extension / on rect(A' 1:s ) such that 

/(¥ = [ d u f( c -u'-y u ) dy u ( 17 ) 

for some point c £ rect(A’ 1:s ). 


7 Scrambled net variance for smooth functions 


Here we prove that the variance of averages over scrambled geometric nets is 
0(n~ 1 ~ 2 ' d log(n) s_1 ), under smoothness and sphericity conditions. The proof 


is similar to the one in Owen (2008) for scrambled nets. We begin with notation 


for some Cartesian products of cells. For this section we assume that dj = d is 
a constant dimension for all j £ l:s. 

Let b be the common base for recursive splits X ? - of X^ C for j £ l:s. 
Let k = (k i,..., k s ) and r = (ti,..., t s ) be s-vectors with kj 
Then we write 

®UKT = Xj,(.kj,tj) ]^[ 


£ N and tj £ 


jeu 


j&u 


and B UKT = rect(B UKT ). For j = 1,..., s, let Sj = ((j — 1 )d + 1 ):{jd) and then 
for u C l:s, define 


Su = U Sj. 

jeu 


(18) 


Now let §„ = {T C S u | T D Sj ^ 0, Wj £ u}. These are the subsets of S u that 
contain at least one element of Sj for each j £ u. There are 2 d — 1 non-empty 
subsets of Sj, and so 

|S„| = (2 d -l)H. (19) 

Lemma 4. Suppose that f is a smooth function on the Sobol’ extensible region 
X 1:s C with extension f. Let each X^) have a convergent recursive split 
in base b whose k-level cells have diameter at most Cb~ k ^ d for 1 ^ C < oo. Let 
u C l:s and let k and r be \u\-tuples with components kj £ N and tj £ , 

respectively for j £ u. Letip UK , T1 be the multiresolution basis function ([7| defined 
by the splits of X 1:s . Then 


I if, 


UKT'yll — 


/ 2\l“l 

( 2 -i) 1 


6-¥(i+l)-^ 


sup \d°f u {y) |. 


( 20 ) 


If f £ C ds (X 1:s ) with Whitney extension f, where each X^ is a bounded closed 
set with non-empty interior and a boundary of measure zero, then (20) holds 
regardless of whether X 1:s is Sobol’ extensible. 
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Note: Recall that we assume that d v f u takes the value 0 in places where it 
is not well defined. Alternatively one could use the essential supremum instead 
of the supremum in (20). Later when we use || • ||oo it will denote the essential 
supremum of its argument. 


Proof. The same proof applies to both smoothness assumptions. From the def¬ 
inition we have 


(/ \lpUKT'y) ukt^y) 

— I I fu (®) dX — u 

J X~ u j X u 

= b~ {W+H)/2 / fu(x) n b k > {b\ Jkjt , c , (xj) - W jkjtj ( Xj )) dx u . 
J X u ■ s- „. 


jeu 


By either Lemma[2]or Lemma[3j f u is smooth and we let f u be its extension. 
We know f u (x) = f u {x) for all x £ X 1:s . As the above integral is over X u , we 
can write it as 


6 -(M+M )/2 f f u{x) JJ b k 3 (p Nk . tjC .( Xj ) _ w kjtj (xj)) dx u . (21) 

^ x " jeu 

Now f u is smooth on rect(A“) and depends only on x u . Applying we can 
write, 


Iu{x) = y l 


CS„ I*"’ 1 ”] 


d v f u (z- v :y v ) dy v , 


( 22 ) 


choosing to place the anchor z at the center of B UKr . Note that if v £ §„, then 
there exists an index j £ u such that Sj D v = 0 and then the integral in p2| 
above does not depend on Xj making it orthogonal to bNj kjtjCi ( Xj). 

Also the integrand in (21) is supported only for x u £ B UKT . Putting these 
together we get, 

6 (W + M)/2| </)V , )| 


E 


ves u 


^ E 


sup 


vGS, 


x u £j. 


d v f u {z- v \y v ) dy v b k > (bN jktc (xj) - W jkjtj (xj)) dx, 

jeu 

d v f u (z- v :y v )dy l 


t X u 


]T[ b 1 | bNjkjtjCjixj) Wjkjtj i x j) | da:,, 


jeu 






sup 

x u £M uk 


d v f u {z- v :y v )dy l 
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Now since d v f u is bounded we can write, 



d v f u (z- v :y v ) dy v 


< vol(rect[z„,®„]) sup \d v f u {y)\. 

3/6 B„„ t 


Because 2 : G B UKT we have 

vol(rect[z„,a:„]) = JJ \z t - x t \ < 

< (C'd 1/2 ) |v| 6 _l ' s|/<i . 

The last inequality follows because v fl Si,- > 1 for all j G u and also uses 
equation (11) on the diameter of a bounding box. Finally, putting it all together, 
we get 

|</,VWr>|< ¥ Y, <?HdH/ 2 sup \d v f u {y)\. □ 


v£$u 


3/6B„ 


The factor d^l/ 2 in the bound can be as large as d s / 2 in applications, which 
may be quite large. It arises as a |u|-fold product of ratios diam(rect(-))/diam(-) 
for cells. For rectangular cells that product is 1. Similarly for cells that are ‘axis 
parallel’ right-angle triangles, the product is again 1. 

Lemma 5. Let h u (z) = max.„ e § u \d v f u (z)\ for z G M UKT . Under the conditions 
of Lemma [^J 


°L< 


If 


i“i 


b~ 2 M /d \\h. 


2 

u 1100 


where C = d 1 > 2 {2 d - \)C d . 

Proof. The supports of ifunr-y and ip UK , T 'y are disjoint unless r = t'. Therefore 

kL(*) = E E</> ^PuKT^y) (f 1 'IpUKT'f , )'4^UKT / y UKT'f' ) • 

t\u 7,7 '!u 

Now 


= 


lx 1 - 


"L( x ) dx 


= E E (/ 5 '^UKT'y) ( f 5 *0 UKT'y' } ^ ^ii/tT7 (*^)V^ , ukt7 / (*^) d.CC 

r|it,« 7 , 7 / |u ^ 

= ^ ^ ' (/) ' t PuKT'y){fl 4’uK.T'l') H(lc,=c' — b 




t\u,k 7,7' |li 

<r 2 - 2 


j 6« 


s«P r/„(y)|) £ n<T«; - 6- 1 )- 

t|u,k, VUGS™ 2/EB ukt / 7,7'! uj£u 
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Some algebra shows that ^ 7 , 7 'l ulljeui 1 ^^ ~b x ) = (2 - 2/6) |u| . The 
supremum above is at most ||/i u ||oo. From equation (19), we have |S M | = ( 2 d — 
and also ^ C d l“l for v £ There are indices r in the sum given 
u and n. From these considerations, 

a 2 UK < (2 - ^ 3|u| 6-2|«l/d|| /lu ||2 o ( 2 rf _ 


€ 


2\ 3 


^( 2 - ) 


1 M 


6 - 2|K|/d ll/i M |lL. 


□ 


Theorem 7. Let u 1 ,..., 6e t/ie points of a randomized (t, m, s)-net in base b. 
Let Xi = (j)(ui ) £ X 1:s for i = 1,..., n where <f> is the componentwise application 
of the transformation from convergent recursive splits in base b. Suppose as 
n —> 00 with t fixed, that all the gain coefficients of the net satisfy T UK < G < 00 . 
Then for a smooth f on X 1:s , 


Var(/t) = O 


f (log n 


|S —1 


\ ^1+2 /d 

If f £ C ds (X 1:s ) where each X^ is a bounded closed set with non-empty interior 
and a boundary of measure zero, then (20) holds regardless of whether X l:s is 
Sobol’ extensible. 


Proof. We know from (10) that 

Var(A) < y E E 

|it|>0 |«i|>(ra— t— |n|)+ 


<- E 

n 2 —' 

| u |>0 

G 


C 2 




M 


W|| 2 oo E 

\K,\>(m—t—\u\) + 


5—2|k| /d 


^E E 


5—2|k |/d 


|u|>0 |ft|>(m— t— |it|)_(_ 


where 


G = G 


2 - 


n M 


max I \h u 

| u |>0 


Since we are interested in the limit as m —> 00 , we may suppose that m > s +1. 
For such large m, we have 

y b - 2 \ K \/d = y' b - 2 r/d( r + M - 1 

\ \u\ — 1 

|/‘t|>(ra— t— |u|) + r=m-t—\u\-\-l 


where the binomial coefficient is the number of |-u|-vectors k of nonnegative 
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integers that sum to r. Making the substitution s = r — m +1 + \u\, 


b~ 2 l K l /d = &C“™+*+M)2/d b~ 2s/d f s + rn ~ t ^ 1 

|«;|>(m—i —|u|)+ s=l ^ ^ 


h (t+\u\)2/d 00 

< b ._ b- 2 */ d {s + m — t— 1)":-' 


n 2 / d (\u\ — 1)! 


S=1 


^i+\u\)2/d 


n 2 / d ([u\ — 


d 00 \ u \~ 1 /I I 

'' s =1 j=0 \ J S 


^ (m-t-l)M---j “ , 

^ jl(M-i-j)! ii 


? 2 /d 


i=o 


j.(t+|«|)2/d _ 

< ^2 7 ^ TOl “ hl i u iE r2s/dslu| “ 1 

S = 1 


Note by the ratio test it is easy to see that b 2s / d sl“l 1 converges. Also 

as to < log h (n) and |u| < s we get 

£ <-- 2M/d =o (i 1 ^ 

\K,\>(m—t—\u\)+ 

Plugging this back into the bound for the variance we get the desired result. □ 


8 Discussion 

Our integration of smooth functions over an s-fold product of d-dimensional 
spaces has root mean squared error (RMSE) of 0(n _1 / 2_1 / d (log(n))^ s_1 ^ 2 ). 
Plain QMC might map [0, l] sd to R. If the composition of the integrand with 
such a mapping is in BVHK, then QMC attains an error rate of 0(n -1 log(n) sd_1 ). 
Our mapping then has the advantage for d = 1 and 2. When the composition is 
not in BVHK then QMC need not even converge to the right integral estimate. 
Then scrambled nets provide much needed assurance as well as error estimates. 

When the composed integrand is smooth, then scrambled nets applied di¬ 
rectly to [0, l] sd would have an RMSE of 0(n~ 3 / 2 log^.)^" 1 )/ 2 ). That is a 
better asymptotic rate than we attain here, and it might really be descriptive of 
finite sample sizes even for very large sd , if the composite integrand were of low 


effective dimension (Caflisch et al. 1997). If however, the composed integrand 


is in L 2 but is not smooth, then scrambled nets applied in sd dimensions would 
have an RMSE of o(n -1 / 2 ) but not necessarily better than that. Our proposal 
is then materially better for small d. 

The composed integrand that we actually use is not smooth on [0, l] s . It 
generally has discontinuities at all 6-adic fractions t/b k for any of the components 
of u. For example in the four-fold split of Figure [1] an e change in u can move 
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a point from the top triangle to the right hand triangle. These are however 
axis-aligned discontinuities. Wang and Sloan (2011) call these QMC-friendly 


discontinuities. They don’t induce infinite variation. 

We have used nested uniform scrambles. The same results apply to other 
scrambles, notably the linear scrambles of Matousek (1998). Those scrambles 
are less space-demanding than nested uniform scrambles. A central limit theo¬ 
rem applies to averages over nested uniform scrambles (Loh 2003), but has not 
been shown for linear scrambles. Hong et al. (2003) find that nested uniform 
scrambles have stochastically smaller values of a squared discrepancy measure. 

The splits we used allowed overlaps on sets of measure zero. We could 
also have relaxed <T = U^~Q<T a to vol(Y \ U^ZgT’o) = 0. That could cause the 
deterministic construction to fail to be a weak geometric ( t , m, s)-net but the 
scrambled versions would still be geometric (t, to, s)-nets with probability one. 

Our main result was proved assuming that all dj = d. We can extend it 
to unequal dj by taking d = rnax^-g!^ dj. To make the extension, one can add 
dj — d ‘do nothing’ dimensions to The splits never take place along those 

dimensions, so the cells become cylinder sets and the function does not depend 
on the value of those components. We can make the extent of those do-nothing 
dimensions as small as we like to retain control of the diameter of the splits and 
then apply Theorem [7j 
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